Nature of the Quantum Phase Transition in Quantum Compass Model 



Department of Physics, 



Han-Dong Chen 
University of Illinois at Urbana- Champaign, 



Urbana, IL 61801 



O 
O 

(N 
O 

Q 



c 

o 
o 



> ■ 
in ■ 
o\ : 

(N ■ 

vo : 
o . 
-i— > ■ 

a ■ 

s; 

T3 ■ 

c : 
o . 
o ■ 



X 



Chen Fang and Jiangping Hu 

Department of Physics, Purdue University, West Lafayette, IN 47907 

Hong Yao 

Department of Physics, Stanford University, Stanford, CA 94305 

In this work, we show that the quantum compass model on an square lattice can be mapped to a 
fermionic model with local density interaction. We introduce a mean-field approximation where the 
most important fluctuations, those perpendicular to the ordering direction, are taken into account 
exactly. It is found that the quantum phase transition point at J x — J z marks a first order phase 
transition. We also show that the mean field result is robust against the remaining fluctuation 
corrections up to the second order. 

PACS numbers: 



INTRODUCTION 



Quantum compass model^ has recently attracted great 
interests^^. It was originally proposed as a simplified 
model to describe some Mott insulators with orbit degen- 
eracy described by a pseudospin. In particular, compass 
model describes a system where the anisotropy of the spin 
coupling is related to the orientation of bonds. More re- 
cently, it was also proposed as a realistic model to gener- 
ate protected qubits^. It has been argued that the eigen- 
states of the quantum compass model are two and only 
two fold degenerate^. The two-fold degenerate ground 
state can be an implementation of a protected qubit if it 
is separated from the low energy excitations by a finite 
gap. Symmetry will protect weak noise destroying the de- 
generacy. Based on the protected qubits, which has high 
quality factor, a scalable and error-free schema of quan- 
tum computation can be designed^. The two-fold de- 
generate state has been shown to be gapped from low en- 
ergy excitations based on the results calculated in a small 
size system^. However, the results can not be extended 
to large size systems. The results from both spin wave 
study and exact diagonalization have suggested that the 
system develops a spontaneously symmetry-broken state 
in the thermodynamic limit 4 . In the spin-wave study, 
the Hamiltonian is expanded around the uniform classi- 
cal ground state up to the first order of 1/S. It is found 
that there is a directional ordering of the ground state 4 .. 
From finite size diagonalization of samples with size up 
to 5 x 5, it is shown that on clusters of dimension L x L, 
the low-energy spectrum consists of 2 L states which col- 
lapse onto each other exponentially fast with L. At the 
symmetric point of J x — J z , 2 x 2 L states collapse expo- 
nentially fast with L onto the ground state. From both 
of the spin-wave analysis and exact diagonalization, a 
first order phase transition at the symmetric point seems 
most favorable. However, the spin- wave analysis ignores 
very important fluctuations while the exact diagonaliza- 



tion is limited by small sample size. In this work, we 
first show that the spin-i quantum compass model on 
a square lattice can be exactly mapped to a fermionic 
model with local density interactions. Normally, after 
performing a Jordan- Wigner transformation, one expects 
a non-local gauge interaction between fermions in the 
fermionic Hamiltonian 8 . However, due to the special 
structure of both the spin interactions and the lattice, we 
show that the gauge interaction for the compass model is 
absent, which allows us to apply conventional approxima- 
tion techniques developed for electron systems to analyze 
the original spin model. Our approximation method au- 
tomatically takes into account the most important fluc- 
tuations, those perpendicular to the ordering direction. 
The remaining fluctuations can also be studied in the 
perturbative approach. It is shown that our conclusion 
is robust against the perturbative corrections. Our re- 
sults support that a first order phase transition happens 
at J x — J z between two different states with spin order- 
ing along either x or z directions. 



II. COMPASS MODEL AND FERMIONIZATION 

The compass model appears rather simple on the first 
look, 

H = — Jx J]] £>i,jSf+l,j ~ Jz J]] S* t jSfj + i, (1) 
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where is a two dimensional coordinate. The sign of 
J x and J z is not important. We can always introduce a 
transformation to bring them into J Xl J z > 0. We assume 
temperature T = throughout this work and set fi — 1 
for simplicity. 

This model has some interesting symmetries, which 
have been discussed in details in literature^. We first 
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recall two types of symmetry generators of this model 3 
Pi=n 2S &> Qj=U 2S ir ( 2 ) 



It has been shown that the symmetry of this model 
leads to a one-dimension type of behavior and direc- 
tional ordering-' 3 ' 5 . It is also shown 5 - that this model is 
dual to recently studied models of p + ip superconduct- 
ing array a 9 ' 10 , which also show the effect of dimensional 
reduction^. 

With the help of the symmetry operators defined in 
Eq.([2]), we find that it is possible to fermionize the full 
compass model. The resulted fermionic model has BCS- 
type pairing along one direction and nearest neighbor- 
ing interactions along the other direction. The Jordan- 
Winger transformation can be defined as 




(3) 



(4) 



where c i j annihilates a fermion at site (i, j). The original 
compass model Eq.([T]) is transformed into 



+ lf ( C ij ~ C L) ( C '+1j + 4+l,j) 



(5) 



up to a constant. Because the interaction along z di- 
rection is quartic, the string-type interaction or gauge 
interaction disappears. 



III. PHASE DIAGRAM 
A. General considerations 

The symmetries defined in Eq.Q also have pro- 
found implications on the possible ordering of this 
model. Based on the fact that [P^P^] = [Qj,Qj>] = 
[Pi, Q3Q3'} = [PiPv,Qj] = and [P, Qj] ^ 0, it is found 
that each eigenstate of this model is two and only two fold 
degenerate^. This leads to the observation that any finite 
system cannot have an ordering that is characterized by 
some finite mean-field expectation value of either or both 
{Sfj) and (Sfj). For instance, let us assume the ground 
state I il) of a L x L system has finite (Sfj). One can 
thus apply Pi onto |0) and obtain 2 L degenerate ground 
state, which is of course inconsistent with the above result 
of double-degeneracy. However, ordering is still possible 
when the system goes to the thermodynamic limit where 
spontaneous symmetry breaking happens. The excita- 
tion gap that separates the true ground state and other 
2 L_1 low energy excitation states collapses exponentially 
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FIG. 1: Four possible scenarios of the evolution of spin or- 
dering, (a) Crossover scenario, (b) Two second order phase 
transitions with coexisting state in between, (c) Two second 
order phase transitions with disordered state in between, (d) 
First order phase transition scenario. 



as the system size goes to infinity 4 . In this case, the 
spontaneously broken symmetries are the Zi symmetries 
of the one-dimensional Ising chain along the ordering di- 
rection. Let us consider a system L x x L z and let both 
L x and L z increase to infinity. If the ordering is along 
the z direction, one should observe 2 Lx_1 states collapse 
onto the true ground state and the gap vanishes exponen- 
tially, namely e~ Lz ^ L ° with some length scale Lq. Since 
the spontaneously broken symmetries are a large num- 
ber of copies of discrete symmetry Z2, there is no cor- 
responding Goldstein mode. Furthermore, one can also 
effectively view this spontaneously symmetry breaking 
as the breaking of a set of local Z2 symmetries in a one- 
dimensional problem, which is along x if the ordering is 
along z. This does not violates the Elitzur's theorem 1 ^ 
since in this case the effective local Z2 symmetry is real- 
ized by a infinitely long Ising chain and thus the energy 
barrier to restore the symmetry is still infinite. 

We shall now consider the phase diagram in the ther- 
modynamic limit. We start from two extreme cases which 
are trivial to solve. We then study how the system 
evolves from one extreme to the other. It turns out that 
by studying this question it is possible to obtain the phase 
diagram as we shall show below. 

Let us first start from the simplest limit J x = 0, where 
the model reduces to decoupled Ising chains. The ground 
state is magnetically ordered such that (S z ) = \ up to an 
overall flip of all spins of each Ising chain. This ordered 
state thus has 2 L degeneracy for a L x L system. Simi- 
larly, at the limit of J z — 0, the spins are ordered along 
x direction. We are thus interested in how the spins ro- 
tate from z direction to x direction as the ratio J z /J x 
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changes from oo to and what is the possible ground 
state at J x = J z . In very general, there are four possible 
scenarios. The first possibility is that there is no phase 
transition. The spins rotate continuously from z to x 
through a crossover. The second scenario corresponds to 
two second order phase transitions with coexisting of or- 
dering along both z and x directions. The third possibil- 
ity is two second order phase transitions with disordered 
state in between. The fourth one is that the spins sud- 
denly change from x direction to z direction, i.e., a first 
order phase transition happens at J z = J x . 



B. Mean- field phase diagram 

To find out which of the four possible scenarios hap- 
pens, we start with limit J x — 0, where the system is 
ordered along z direction, and increase J x (or decrease 
the ratio J = J z /J x ). As J x is increased, the flipping 
of two spins on adjacent chains is introduced. This pro- 
cess tends to reduce the ordered moment along the z- 
direction. On the other hand, the ordering along z direc- 
tion acts as an effective transverse field and suppresses 
the ordering along x direction of the spins on adjacent 
chains. This effect is described by a ID Ising chain with 
a transverse field, which has a critical transverse field 12 . 
When J x is much smaller than J z , the ordering of spins 
along x direction is fully suppressed by the adjacent spins 
ordered along z direction. We then ask at what value of 
J x the spins start to order along the x direction. For 
the crossover case, (S x ) has a finite value as long as J x is 
not zero. For the second scenario, two second order phase 
transitions with intermediate coexisting state, (S x ) starts 
to develop at some value of J that is larger than 1. For 
the last two cases, the expectation value of S x is zero for 
all J x < J z , while (S z ) takes a finite value for all J z > J x 
in the first order phase transition scenario. To determine 
the behavior of (S x ) and (S z ), we self-consistently study 
the mean-field decoupled Hamiltonian 

Hmf/Jx = - J~] SfjS? +1 ^ - B ef f (ijSlj , (6) 



B eff (i) = J + (S* j+2 )} 



(7) 



The mean-field decoupling is justified, since we start from 
the situation where S z is ordered and want to know when 
S x starts to develop a finite expectation value. This mean 
field decoupled Hamiltonian describes the most impor- 
tant fluctuation in the state with spins ordered along z 
direction. The ordering along x direction is suppressed 
by the ordering along z direction, which behaves as an 
effective transverse field for S x . On the other hand, the 
feedback effect on S z due to S x fluctuations is also in- 
cluded in this approach through the self-consistent con- 
dition 

(s? d ) = (m,M, (8) 

where \ip) is the ground state of the mean-field Hamilto- 
nian Hmf- To minimize the ground state energy, it is 



favorable to have (S"j) j-independent, i.e., the moment 
on the same chain is uniform such that the J z terms are 
minimized. The mean-field Hamiltonian thus describes a 
series of decoupled Ising chains with an effective trans- 
verse magnetic field B e ff = 2J(S'f J ) that is determined 
self-consistently. 

To minimize the coupling energy along z-direction, 
one would expect \B e ff(i)\ = B e ff is i-independent. If 
B e ff(i) = —B e ff, we can rotate the spins on i-th row by 
an angle 7r around the x-axis. Without losing generality, 
we can thus consider the case where B e ff(i) — B e ff. 

In the Fermionic language, the mean-field decoupling 
introduced in Eq.© corresponds to the decoupling of 
the interaction term J z n i . +1 - 
chemical potential term pm^ with 



J z n, „• into an effective 



H = 2J{Uij) - J = B e ff. 



(9) 



This Hamiltonian can be diagonalized using Bogoliubov 
transformation 

Ik = u kc k + v k d k 7_ fc = v k c\ - u k c_ k , (10) 

where c k is Fourier transformation of c(i). The u k and 
v k can be solved to yield 




(cosfc)/2 — \i 



E k 



^ (cosfc)/2 — fi 
E k 



(11) 
(12) 



with E k = ^Jl/A — pcosk + /i 2 . It follows that the ex- 
pectation value of S z is 



2tt 



2tt 



2B P f r — cos k 



8E k / J x 



The self-consistent equation thus reads 



B eff ~ 



dk 



cos k — 2B, 



ff 



1 - AB ef f cos k + AB 



(13) 



(14) 



eff 



From the results shown in Fig[21 we know there are 
nontrivial solutions to the self-consistent equation (fT4]) 
for J > 0.7446. The jump happens at i? e // ~ 0.65 > 
B c = 1/2. Below B c , which is the critical point of the 
Ising chain in a transverse magnetic field, S x starts to 
develop a nonzero expectation value. We notice two im- 
portant results that can help us determine the phase 
diagram. Firstly, the critical point B c happens at a J 
smaller than the symmetry point J = 1, below which 
ordering along z direction is no longer a reasonable as- 
sumption. For J < 1, the ordering along x direction is 
instead more favorable. This implies that (S x ) cannot 
develop a finite value for J x < J z . The first two sce- 
narios, namely crossover scenario and two second-order 
phase transitions scenario, are thus excluded. Secondly, 
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Hamiltonian Eq.© 




FIG. 2: The solution to self-consistent equation (|14|l . There 
is no non-trivial solution for J < 0.7446. (a) The effective 
transverse field B e // vs J. (b) The ordered moment along z 
direction (S z ) vs J. 



the self-consistent equation has non-trivial solutions for 
J > 0.7446, which means spins order along z direction 
for J z > J x . Combining these two points, we conclude 
that there is a first order phase transition at J = 1. At 
J > 1, the spins are ordered along z direction while they 
are along x direction for J < 1. At J = 1, the spins are 
either along i or z direction since our argument above 
shows that the ordering at J = 1 along x (z) direction 
induces a strong enough transverse field to suppress the 
ordering along z (x) direction. This conclusion is consis- 
tent with the spin- wave analysis*. In spin wave analysis, 
fluctuations of both directions are taken into account par- 
tially. In our approach, the most important fluctuations, 
namely the ordering direction of weaker bond, is solved 
exactly. 



Hi, 



= Jznijriij +1 + 2J z (n) 2J n 
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(15) 



In terms of Feynman diagrams, there are hence two 
kinds of vertices: two-particle interaction vertex and ex- 
ternal field interaction. Let us first work out the Feyn- 
man rules for the two-particle interactions. We observe 
that the interactions are between particles on nearest 
neighboring lines. The vertex of two-particle interaction 
thus consists of three parts, the coupling strength, the 
contribution from j-th line and the contribution from the 
neighboring (j + l)-th or (J — l)-th line. There are four 
possible types of contributions, which are given in Fig[3] 
with 



tanWfe = 



— sin A; 



2J Z /Jx + cos k 



(16) 



<: 
k. 



:cos(e k2 +e k1 )/2 = B + .(k 1 ,k 2 ) 
= -cos(9 k2+ e k1 )/2 = B + (k 1 ,k 2 ) 



dsin(e k2 -e ki )/2 = B ++ (k i ,k 2 ) 



:isin(ft 1 -ej/2 = B..(k„k; 



FIG. 3: Four possible contributions of j-th line to the vertex. 

The first line corresponds to (7fc 2 7fc 2 )(7fc 1 7fc 1 )- In this 
section, we will use (A) to denote the expectation value 
of operator A in the un-perturbed mean-field ground 
state. This term can be obtained by contraction from 
either <7fc 2 (7l 2 7/^)7^) or ( 7fe ( 7fci 7^ ) . The coher- 
ence factor of (7^ 7fc ) in the two-particle interaction is 

cos cos ^ while the one of {"i kl l\ 2 ) is sin sin 
Taking into account of the sign difference between these 

two contractions, we find (fci,^) = cos kl 2 k2 . 

Similarly, we get the results for other three possibilities. 



C. Beyond mean-field approximation 

Although the mean-field approximation has included 
the most important fluctuations, there is still concerns 
about the effect of the ignored fluctuations. To fur- 
ther study the effect of fluctuations, we calculate the en- 
ergy gap by perturbation. The perturbation comes from 
the interaction terms that are ignored in the mean-field 




= - JzB , (k 1 ,k 2 )B lt1 (k 3 ,k 4 ) 



FIG. 4: The vertex of two-particle interaction. 
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FIG. 7: First order correction to the self-energy. 



FIG. 5: The vertex for the external field term. 
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FIG. 6: The Feynman rules for bubble and propagator of 
particle. 



The resulted vertex for two-particle interaction is thus 
—J z Bi{ki, k2)B : ' +1 (k%, ki) with corresponding Feynman 
diagram sketched in FiglU Here, we have not shown 
the arrows associated with ki,k2,ki + q, fc 2 — Q- These 
arrows determine the corresponding subscripts of _B- J+1 
and B J . The corresponding momentum conservation law 
for the vertex is S(±ki ± k 2 ± k^ ± k^) where the signs 
of ki(ks) and ^2(^4) are the first and second subscripts 
of B^B^ 1 ). For instance, the momentum conservation 

for B 3 + _(ki, fc 2 )i?^."t 1 (fc3, ki) is 5(k\ +k 3 — k 2 ~ k^). Simi- 
larly, we can obtain the vertex for the external field term 
as illustrated in FigJS] 

We are now at position to calculate the Feynman rule 
for bubbles. The clockwise bubble corresponds to (j^l) 
with coherence factor sin 2 ^ while the counterclockwise 

one is (7^7 fc ) = with coherence factor cos 2 ^. The 
Feynman rules for bubble and propagator of particle is 
shown in Fig [6] A noticeable difference between this prop- 
agator and the one in the familiar case of electrons under 
coulomb repulsion is that we only have the plus sign be- 
fore the infinitesimal imaginary part. This is because in 
the current system, we don't have a Fermi sea as ground 
state and therefore do not have a hole propagator. This 
property of system largely reduces the number of dia- 
grams we need to calculate. Finally, every intermedi- 
ate 4-momentum should be integrated with a factor of 
H)/(2tt). 

With all the Feynman rules in hands, we are now ready 
to do perturbative calculations. The first order correc- 
tion to the self-energy is given in Fig(7] The first term is 
the contribution from the interaction with nearest neigh- 
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FIG. 8: Second order correction to the self-energy. 



bor lines. The factor 2 comes from the fact that there are 
two nearest neighbor lines. The second term is from the 
"external field" term in Eq. (|15p . The first order pertur- 
bation vanishes in the mean- field decoupling of Eq.([S]). 
The leading order correction is thus second order per- 
turbation given by the connected amputated graphs as 
shown in Fig|5J 

There are five pieces of nonzero second order correction 
to the single particle propagator. We detail the calcula- 
tion of them separately. The first two diagrams denote 
the self-energy process involving two particle-particle in- 
teractions. The overall factor of two stems from the fact 
that a particle on the, say, jth chain can either interact 
with the (j — l)th or the (j + l)th chain. The relative 
minus in the bracket is there because we invert three (an 
odd number) of the propagators in the second diagram. 
Direct application of the Feynman rules we derived yields 



h 



h = 



2tt COS' 



i,J z 



(^4^) sin 2 ( ii^i^z£i2. \ 

W — £fei — £/c 2 — e fe-fci-fc 2 



dk\dk2 
(17) 
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FIG. 9: (a) Energy gap of the ground state with second order 
perturbation, (b) The ratio between second order correction 
and the energy gap itself. 

• 2 (Ok— 9k , \ • 2 f 8-k-k-, -k 2 -0 k2 \ 
r2n Sin I 2"^- 1 Sill I ^y- 2 2- 

x / i '- i J -dkxdk 2 , 

Jo W + Efej + £fe 2 + £fe_fe 1 _fc 2 

(18) 

where I\ and -Z2 are the contributions from the first and 
second pieces of Fig [8] respectively. The intermediate fre- 
quencies oji and UJ2 have been integrated. 

The third and the fourth diagrams denote self-energy 
processes involving one particle-particle interaction and 
one external field perturbation. The factor of 2 stems 
from the same argument and the overall negative sign 
can be obtained by working out the contractions. Their 
analytical expressions I3 and I4 read 

h=h = -2J z (n) j ^ — dk , (19) 

where again, the intermediate frequency has been inte- 
grated. 

The last piece of the expression denotes process that 
involves two external field perturbation. It is easy to 
obtain the analytical expression as 

h = -i{2J z {n)f S ^ k - (20) 
lu + e k 

In this way, we have exhausted all possible self-energy 
processes to the second order of J z . The resulted sum of 
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FIG. 10: (a) The ordered moment along z direction versus 
Jz/Jx- (b) The effective field -B e // = 2J Z S Z after second- 
order correction is included. Solid lines are the results with 
second order perturbation included while dash lines are mean- 
field results. 



the diagrams is given by — iE(w, k) = 2I\ — 2I 2 — 2/3 — 
2Z4-I-I5, where £(<x>, k) is the correction to the propagator 
of a quasi-particle carrying momentum k. The energy of a 
particle can be defined as the singularity of its propagator 
and readily obtained by solving the following equation 



E(k)-e k -H(E(k),k) = 0. (21) 



The energy gap is defined as the lowest energy of excita- 
tion, which is numerically calculated and shown in FigJ5] 
The second order correction amounts to about 30% of 
the corrected energy gap itself at the symmetric point 
Jx — Jz- In FigfTOl we also plot the ordered moment S z 
as a function of J z /J x - Again, the moment is reduced 
by a small second-order correction. At the symmetric 
point, the effective field is still larger than the critical 
value B c = 1/2 below which an ordered moment along 
x direction becomes possible. We thus conclude that the 
mean-field results are robust against the second order 
perturbation. 
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IV. CONCLUSION 

In conclusion, we show that the compass model can be 
mapped to a fermion model with local density interac- 
tion. Through a self-consistent solution of the model, we 
argue that there is a first order phase transition at the 
symmetric point J x — J z . This conclusion is consistent 
with spin- wave analysis and recent numeric computations 
of the spectrum^. In our approach, the most important 
fluctuation, namely the fluctuations along the weak cou- 
pling direction, is taken into account exactly. The fluc- 
tuations along the ordering direction is considered up to 



second order perturbation. It is shown that the result 
of our mean-field approximation is robust against such 
corrections. 
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